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Abstract 

We present calculations of the NN — > NA cross section in nuclear matter within the one- 
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nuclear interactions and Dirac effective masses of the baryons due to coupling with the scalar a 
field. Introducing the Dirac effective masses leads to an in-medium reduction of the cross section. 
The experimental data on pion multiplicities from the collisions of Ca+Ca, Ru+Ru and Au+Au 
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I. INTRODUCTION 



While transport calculations have, in general, been very successful in describing particle 
production in heavy-ion collisions [l|, , one of the remaining 



open problems in the transport-theoretical description of heavy-ion collisions at SIS energies 
(~ 1 A GeV) is an overprediction of the pion multiplicity At the 

beam energies of a few A GeV the dominating mechanism of the pion production is the 
excitation of the A(1232) resonance in a nucleon-nucleon (NN) collision followed by its decay: 
NN — > NA, A — > Ntt. The pion multiplicity, therefore, depends crucially on the value of 
the in-medium NN — > NA cross section which - by detailed balance - also determines the 
pl on reabsorptio, In Ref. Q a phe—logical dependent —ion facto, 
of the vacuum NN — > NA cross section has been proposed. The BUU calculations with 
the density dependent NN <-> NA cross sections are in a reasonable agreement with the 
experimental data on pion production |l3| |. 

The purpose of the present work is to get rid of the ad-hoc description of Ref. jl^ and 
to calculate the in-medium NN <-» iVA cross sections on the basis of the one-pion exchange 
model (OPEM) with medium modifications. The calculated cross sections are used then in 
the BUU model to study the pion production in heavy-ion collisions at SIS energies. 

It has been shown in Ref. 3] within the relativistic Dirac-Brueckner approach that the 
NN — > NA cross section is reduced at high densities. Later studies within the OPEM 
including pion collectivity and vertex corrections Q, Q] have shown, however, that the 
NN — > NA cross section grows with nuclear density. An important contribution which 



was missed in Refs. 



:ction gr 

HQ 



171 ] is given by the Dirac effective nucleon and A resonance 



masses. We will show, that the introduction of the Dirac effective nucleon and A masses in 
the OPEM leads to an in-medium reduction of the NN <-> iVA cross sections in agreement 
with Ref. 

The structure of the paper is as follows: In Sect. II we describe the OPEM with medium 
modifications including pion collectivity, vertex renormalization by the contact interactions 
and effective masses for the nucleon and A resonance. The energy dependence of the total 
pp — > nA ++ cross section as well as the angular and A resonance mass dependencies of the 
differential pp — > nA ++ cross section are studied at finite nuclear matter densities. Sect. Ill 
explains how the in-medium modified cross sections are implemented into the BUU code and 



2 



describes the results of the BUU calculations of pion production from heavy-ion collisions at 
0.5-=-2 A GeV. In sect. IV the summary of our results and of remaining problems is given. 



II. THE MODEL 

For the calculation of the iViV — > NA cross section in nuclear matter we apply the 
nonrelativistic OPEM similar to Ref. jlfij . however, with readjusted parameters. We have 
chosen the nonrelativistic version since it incorporates the contact nuclear interactions in a 
natural way, which is necessary in the in-medium calculations. 

In the nonrelativistic reduction (c.f. Ref. the irNN and ttNA interaction La- 

grangians are: 

f 

£kNN = WGaT 1pV a 7Z , (1) 



ArVA = — 4> A S a Tip\7 a n + h.c. , (2) 
m 



where ip, ip A and 7r are the nucleon, A resonance and pion field respectively. / = 1.008 and 
/a = 2.202 are the coupling constants 19]. a a and r are the spin and isospin Pauli matrices. 
S a and T are the spin and isospin transition (1/2 — > 3/2) operators defined according to 
Ref. IJ. 

The nuclear spin-isospin short range correlations are known to be important for the in- 
medium iViV — > iVA cross section We will introduce them via the following Lagrangian: 

C-SRC = -W-ZSnnW VaT l/j)^ (T a T if)) 



m 2 



P 

J A J q rp„/A/ a/l t 



r ^ AA (^S , a TV)(^«T^)+h.c. 



_2ml ' 

9nni 9na an< i 9aa are ^e Landau-Migdal parameters. 



(3) 



The values of the Landau-Migdal parameters are not fixed unambiguously in the liter- 
ature. Within a simple universalit y a ssumption g' NN = g' NA = g' AA = g' BW , which is the 
Backmann-Weise choice (c.f. Ref. 21] and Refs. therein), one gets g^iy — 0.7 ± 0.1 from 
the best description of the unnatural parity isovector states in 4 He, 16 and 40 Ca. However, 



the same calculations within the Migdal model |22| assumption g' NA = g AA = produce 



9nn = 0-9 1- The description of the quenching of the Gamow- Teller matrix elements 



requires g' AA = 0.6 4- 0.7 (assuming g' NA = g AA ) [23[. From the real part of the pion optical 
potential in 7r-atoms one gets g' NA = 0.2 and g' AA = 0.8 |22j. The pion induced two-proton 
emission is the best described with g' NA = 0.25 4- 0.35 |24|. We will adopt the two sets of 
the Landau-Migdal parameters from Ref. [25j (see Table 1). 

As a first step, we calculate the NN — > NA cross section in vacuum. The in-medium 
corrections are introduced in a second step. 

A. Vacuum cross section 



The differential N1N2 — » A3A4 cross section can be written as: 

da = (2n) 5^(p 1+ p 2 -p 3 -p,)\T\i ^ {27[?2£ ^ {M A )dM A , (4) 



where |T| 2 is the square of the absolute value of the reaction amplitude T in the normalization 



of Ref. 



26| averaged over spins of incoming particles 1,2 and summed over spins of outgoing 



particles 3,4. The flux factor is : / = y (P1P2) 2 — m %- The spectral function of the A 
resonance 

a (M 2, 1 M A r A (M A ) 

where m A = 1.232 GeV is the A pole mass, satisfies the normalization condition : 

00 

J dM 2 A A A (M 2 A ) = 1 . (6) 

(m N +m n ) 2 

The mass dependent total A width Ta(Ma) is parameterized according to Ref. : 

r A (M A ) =r° ( g( M A,rajv,wyU 3 m A fil + g 2 (m A} m N} m n ) 
A A A \q{m A ,m N ,m n ) J M A ffi + q 2 (M A ,m N ,m n ) 

where T° A = 0.118 GeV is the width at the pole mass, /3 = 0.2 GeV/c is the cut-off 
parameter. In Eq.((7j) and below 



. . Urn 2 + m\ - ml) 2 , , oN 
q(m,m 1 ,m 2 ) = y m( (8) 

is the center-of-mass (cm.) momentum of outgoing particles with masses mi and 1712 from 
the decay of a particle with mass m. In the cm. frame of colliding nucleons one gets after 



standard transformations : 

<fr (2m JV ) 3 2M A -— g(y/s,M^,m N ) . 2) , . 

where s = (pi +]?2) 2 - 

The amplitude T is evaluated from the graphs shown in Fig. ^ In order to preserve the 
symmetry of the cross section with respect to Q c .m. — 90° we have included the contact 
interaction both in the direct and in the exchange diagrams. The same was done in Refs. 



15. 



12811 whereas in Ref. 



24j the contact interaction has been included in the direct term 
only. For the discussion on this subject we refer the reader to the Appendix of Ref. 
Applying Feynman rules with the Lagrangians (0)-© one gets the following expressions for 
the direct (a) and exchange (b) terms of the amplitude : 

iT a = lx ] A (4:)rf s (k)x(2)]iD {0 \k){x^(3)Tf c (k)x(l)] - iV 4>3;2 ,i{k) , (10) 
iT h = -[xL(4)rf {k')x{l)]iD^\k>)[x\Z)Tf\k')x{2)\ + iV m (h>) . (11) 

Here x an d Xa are the spin-isospin wave functions of the nucleon and delta : 

X — $m t mfi\s j XA = $M t M$\ A s A > (12) 

where m t = ±1/2 (M t = ±3/2, ±1/2) and A = ±1/2 (A A = ±3/2, ±1/2) are the quantum 
numbers of the third isospin and spin component of a nucleon (delta). The quantities m, 
s, M and s A are the corresponding indices of the wave functions. The four-momenta of 
the exchange pion in the direct and exchange channels are k = p\ — p% and k! = p 2 — p%, 
respectively. Below we use the kinematical invariants t = k 2 and u = k' 2 . The operators 



Tf s (k) = -i^-(S • k)Tj , (13) 



m 

Yf%k) = ^(<x-k)r, (14) 
represent the ttN — > A absorption and iV — > ttN decay vertices in vacuum. The quantity 

v mi {k) = ^PiAkiwwii^^^i)] (is) 

7T 

is the contact interaction in vacuum and 



is the free pion propagator. In the ttNN and 7riVA vertices we use the monopole form factor 
F : 

A2 _ rr) 2 

fit) = f F{t), Uit)=f A F(t), F{t) = -^—^ (17) 



3. [uJQ 



with A being the cut-off parameter (see Table 1). Following Refs. |l5l . l20j a factor F 2 is 
included in the contact interaction vertices, since these vertices simulate the exchange by 
some heavy meson. 

The matrix element squared is 

\T\ 2 = \T a + T b \ 2 = \T a \ 2 + \T b \ 2 + T a T* + T a *T 6 . (18) 

We now consider the case of the pp — > nA ++ reaction. The spin-averaged direct term in 
flUJ) is 

2 

/JVA 



\ E W = U^P^) [k^ 2 {k) + 2k 2 D^\k) g ' NA + 3ig' N 

Ai,A 2 ,A 3 ,Aa4 V 7T / 



2 ] • (19) 



The spin-averaged exchange term ^ ^2\i,\ 2 ,x 3 ,x^4 l^~fe| 2 is given by Eq. (|19|) with replacements 
k — ► fc' and t —>■ u. The spin-averaged interference term can be written as follows : 



\ e ew + T:r b ] = - /A(t)/(t)/ t M/M |((k ■ ko 2 + k 2 k' 2 ) J D(°)(fc) J D(°)(^) 

Ai,A 2 ,A 3 ,A A4 2m 7r 16 



(20) 



In numerical calculations we also took into account the cm. and relativistic corrections 
to the amplitudes GD^flDJ (cf. Ref. 3) : 0) We replaced in fjTT^l (S • k) -» (S ■ k 4 ) 
and (S • k') — > (S ■ k^) in the direct and exchange terms respectively. Here k 4 (k^) is the 
momentum of the exchanged pion in the rest frame of A for the direct (exchange) term : 

(ii) The direct amplitude T a is multiplied by the factor (— t/k 2 ) 1 / 2 and the exchange ampli- 
tude T b by the factor (— -zz/k' 2 ) 1 / 2 . This restores approximately the Lorentz invariance after 
summation of the amplitude squared over spin projections. 

After the cm. and relativistic corrections Eqs.(|19 |l . (j2()|l get slightly modified. We do not 
show here explicit formulas: they will be given below for the more general in-medium case 
(see Eas.flA"7l). (lA"8j) in Appendix A). 



In the upper panel of Fig. El we present the total cross section <r pp ^ n A++ as a function 
of the cm. energy. Dashed and solid lines correspond to the case of Set 1 and Set 2, 
respectively; the dotted line represents the relativistic calculation using the model of Ref. 
|19j | . The relativistic and our corrected nonrelativistic calculations are in a good agreement 
at y/s < 3 GeV. At higher cm. energies our calculations produce somewhat larger cross 
section than the relativistic result of Ref. In the lower panel of Fig. El we show the total 
cross section of the process cr pp ^ pnn + vs the cm. energy in comparison to the data from 
Ref. j29|. In this case the dominating channel is still cr pp ^ n &++, A ++ — > pn + , however, also 
other channels contribute: <j pp ^ p a+, A + — > mr Jr , which has a cross section of 1/9 of the one 
for the channel through the A ++ production; channels with higher resonance excitation 
whose contributions are shown by long-dashed and dash-dotted lines for the isospin 1/2 and 
3/2 resonances, respectively; the nonresonant (background) s-wave channel pp — > npn + j]| 
shown by short-dashed line. The sum of all contributions in the case of our pp — > nA ++ cross 
section calculated with the Set 2 and in the case of the cross section from Ref. ^| is shown 
by the solid and dotted lines, respectively. The experimental data are well described in both 
cases. For ^Js < 3 GeV the Set 2 is a little closer to the data, while for higher cm. energies 
the relativistic calculation of Ref. works somewhat better. For the practical applications 
the difference at high cm. energies is not important, since anyway the FRITIOF model is 
applied for the baryon-baryon collisions at ^ > 2.6 GeV in the BUU program (^tJ which 
we use in the description of pion production from heavy-ion collisions. 

Since our major goal is to describe the pion production from heavy-ion collisions at the 
beam energies ~ 1 A GeV, we have also checked the differential cross sections pp — > np7T + 
at piab = 1-66 GeV/c (Ei ab = 0.97 GeV). Fig. El shows the cm. polar angle distribution of 
neutrons. Our calculation is again in a very close agreement with the relativistic treatment 
of Ref. [l^l and describes the data [3^ well. We neglected the contribution of the channel 
pp — > pA + , A + — > 7i + n in comparison to the data: taking into account this channel would 
roughly multiply our results by a factor of 10/9, which would yield an even better agreement 
with the data. The invariant mass distribution of (n + ,p) pairs is shown in Fig. 0] We see 
that both models agree quite well, but the experimental distribution from (^(J is somewhat 
broader. Notice, that the channel pp — > pA + , A + — ► 7i + n missed in our calculations would 
create some broadening of the spectrum, since now the pair (tc + ,p) is not emitted from the 
same resonance. However, the study of this effect is beyond the scope of our work. 
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B. In-medium corrections 



1. Effective masses and in-medium A width 

We assume that nucleons and A's are coupled to the scalar mean field o by the same 
universal coupling constant g a [31]. This produces the Dirac effective masses 

m* N = m N + g a a , m* A = m A + g a a , (23) 

which are smaller than the bare ones by 300-=-400 MeV at the normal nuclear density, as 
one can see from F ig. El We have adopted the two parameter sets of the relativistic mean 
ield model (RMF) |32|: NL1 and NL2, and the Relativistic Hartree Approximation (RHA) 
. RHA takes into account vacuum fluctuations which are missed in RMF. We observe 
from Fig. EJ that the RMF set NL1 produces the steepest decrease of the nucleon effective 
mass with the baryon density, while RHA gives the slowest one. 

Assuming, further, that the coupling constants to the vector mean field us are also the 
same for nucleons and A's 31], we get the kinetic four-momenta 

P*n=Pn~ gj*> , P*a=Pa~ gj*} • (24) 

The effective masses and the kinetic four-momenta substitute the vacuum masses and four- 
momenta in the calculations. Thus, in Eq. (0J) one should replace : 

Pi ^p*, i = 1,2,3,4, (25) 
m N -> m* N , (26) 
M A M* A , (27) 



I ^I* = J(pM) 2 -m**, (28) 



A A (Ml) a* a (m a 2 ) = ^^.Jg^V^) • (29) 

The total in-medium A width is approximated as follows : 

F* l M* \ — T P 1 F° ( <l( M L m *N, m n) Y m* A ffi + q 2 {m* A ,m* N ,m„) 
A[ A) Lsp Po + A {q(m A ,m N}m7I )) M* A f% + q^M A ,m* N ,m n ) ' [M) 

where po — 0.16 fm~ 3 is the nuclear saturation density and T sp is a constant. The first term 
in Eq. (|3T?|) is the A-spreading width 0, [3^ , which describes the in-medium broadening of 
the A resonance due to the absorption and rescattering processes AiV — > NN, AN — > AN 
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and ANN — > NNN. In Ref. |35( the many-body calculations of the A self-energy in nuclear 
matter have been performed which are in a good agreement with the empirical A-spreading 
potential of Ref. [34( 

W sp = (23 ±5) -i(43±5)y (MeV) . (31) 

Thus, the value of the constant T sp consistent with Refs. jsJ, [3| is T sp = — 2lmW sp (p = 
Po) = 80 MeV. In Ref. 3(| another value T sp = 20 MeV has been used in the calculation 
of the pion self-energy within the A-hole model. In Ref. [stJ the A-spreading width has 
been calculated in the relativistic meson-nucleon model including effective masses given by 
the Walecka model. Only the AN — > NN contribution has been taken into account. Vertex 



corrections were neglected in Ref. [37j, which resulted in rather big values of the A-spreading 
width in the calculation without effective mass (T sp = 160 MeV). Including effective mass 
corrections, the A-spreading width of Ref. js^], first, increases with density at p < p in 
agreement with T sp = 80 MeV, but than drops abruptly at higher densities due to the mean 
field effects. We have checked that the variation of the parameter T sp within the range -i-80 
MeV leads to the changes in the in-medium cross section NN — > NA which are comparable 
with the ambiguity due to the choice of the effective mass parameterization (NL1 or NL2). 
On the other hand, we need to keep T sp ^ in order to have numerically stable results in 
the calculations of the A-hole Lindhard function f)44j) . In the numerical results presented 
below we use T sp = 80 MeV. 

The second term in Eq. (j30j) is the A — > Ntt width taking into account the effective masses 
of the nucleon and A. In this term we have omitted the Pauli blocking of the decay nucleon. 
This is justified because in the real heavy-ion collision at the beam energy of 1 A GeV the 
nuclear matter at finite temperature of about 70 MeV is created, which reduces the Pauli 
blocking effect strongly with respect to the zero temperature nuclear matter. 

In the cm. frame of colliding nucleons (i.e. in the frame where P1+P2 = 0) the differential 
cross section in nuclear medium reads as follows (c.f. ©) : 

da^ (2m^) 3 2MA^ ¥ g(y^,MA,m- v ) ^ ( ^ 2) ^ 



dMfdVL 64vr 2 s* 1 1 q(y/^,m* N ,m* N )' 

where s* = (p\ + p^) 2 - Below we will assume that the cm. frame of colliding particles 
coincides with the nuclear matter rest frame, where spatial components of the vector field 
disappear, thus p* = p. 
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2. Pion self-energy 



The in-medium matrix element T med in Eq. (|32|) is given by the graphs shown in Fig. 
where now the in-medium pion propagator 

D{k) = ^- m ;-n W (33) 

must be put on the place of the vacuum propagator D^°'(k) and the in-medium vertex 
corrections due to the contact interactions must be also taken into account. The pion self- 
energy Yl a b(k) is, in general, a matrix with two isospin indices a, b = 1, 2, 3. However, in this 
Section we consider isospin-symmetric nuclear matter, where U a b(k) = 5 a i,Il(k). The delta 
function of the isospin indices will be always dropped below for brevity. The assumption 
of isospin-symmetric nuclear matter results, in particular, in the validity of the vacuum 
relations between various isospin channels (c.f Eqs. (f5T)|) . (j57j) below) for the in-medium cross 
sections. To avoid misunderstanding, we stress, however, that in the BUU calculations of 
Sect. Ill we take into account the isospin assymmetry of colliding nuclei explicitly. 

The function —iU(k) = —iU A (k)—iIl^(k) is shown in Fig.|HK It includes the A-hole and 
nucleon-hole iterated contributions: we denoted them as — iU A (k) and — iIIjv(A;) respectively. 
These contributions satisfy the recurrence relations shown in Fig. |Hb,c which read as : 

/ 



A(k) 

where the 2x2 matrix A(k) is 



ll N (k) = _ k2 XN(k) 
\U A (k))~ \ X A(k) 



(34) 



A(k)=\ 1 + | . (35 ) 
g' N AXA(k) l + g' AA XA(k) 

The susceptibilities Xw(&) an d Xa(^) are : 

XN (k) = -±\+<l>(k) , (36) 
mi 

XA(k) = 16f Q ^f ] (Mk) + M-k)) (37) 
with 4>(k) and 4> A (±k) being, respectively, the nucleon-hole and A-hole Lindhard functions : 

m = i J ^G(P)G(P + k) , (38) 

<j> A (±k) = i J 0- 4 G(p)G A (p±k) . (39) 
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Here G(p) and Ga(pa) are the nucleon and A propagators : 



1 

p° - £*(p) + zO 



G (P) = o I ,n + 2tt» n(p)(y(p° - e*(p)) , (40) 



Ga(pa) = , (41) 

P°A-^A(PA)+^( v /pi)/2 

where 

^p) = ypHKF, ^a(pa) = VpI + (™a) 2 (42) 

are the on-shell nucleon and A energies, n(p) is the nucleon occupation number. The 
nucleon propagator (jjQ)) consists of the vacuum part and the in-medium part (oc n(p)). 
The A propagator (f4*Tj) includes the vacuum part only, since we have neglected the presence 
of A excitations in nuclear matter. Both propagators take into account the effective mass 
corrections. 

I — I 

After the contour integration over p° (c.f. |38|) the nucleon-hole Lindhard function (|38|) 
takes the following form : 

, rA A f d3 P ( Tl(p + k)(l-7l(p)) 7l(p)(l-7l( P + k)) \ 

9{ ' J (2n) 3 [s*(p + k)-k°-e*(p) + i0 e*(p) + k°-e*(p + k) + iOj ' { ' 
We have used the analytic formulas given in Ref. [38] for the calculation of the integral 
in the case of Fermi distribution at zero temperature n(p) = 0(pf ~ |p|)- 

The A- hole Lindhard function (J39|) . after integration over p°, can be expressed as 

M±k) = - ( ^ , (44) 

J (2vr) 3 g*( p ) ± k° - e* A (p ± k) + iT* A ( J(p* ± k) 2 )/2 



where p* = (e*(p),p). In the calculation of the integral ()44j) we have used the expression 



from the Appendix of Ref. 



(see also the Appendix in Ref. j^) assuming that the 



invariant mass of the nucleon and pion in the argument of T* can be approximated by 



(p* ± k) 2 ~ (m* N ) 2 + k 2 ± 2k°^j (m* N ) 2 + -pf, . (45) 

Solving the system (|34j) we get for the pion self-energy the following expression : 

n(Jfe) = U A (k) + U N {k) = -k 2 det _1 (i(A;)) 
x [(l + g'AAXA(k))xN(k) + (I + g' NN XN(k)) X A(k) - 2g' NA x N (k)xA(k)] , (46) 

where 

det(i(A;)) = (1 + 9NNXN(k))(l + S'aaXa(^)) - G4a)V(£0xa(A0 . (47) 



Expression (j4*U|) coincides with the result from Ref. 41 1 
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3. Vertex corrections 



The renormalized ttN — > A absorption vertex operator r| 6s (/c) is shown in Fig. EH- The 
first graph in Fig. [7^ is the 7riV — > A absorption vertex operator in vacuum Tj bs (k) (see 
Eq.()13j)). The second and third graphs in Fig.[7k correspond to the iterated nucleon-hole and 
A-hole corrections which we will denote as T^(k) and ■ ^(k) respectively. The correction 
terms satisfy the recurrence relations shown in Fig. [?b,c which can be written as : 



A(k) 



ft 



W | = _ rabs[k) ( g' NA XN(k) \ 
\r?i(k) ) 3 {g' AA XA(k) ) 



(48) 



After solution of the system (J48|) we obtain the expression for the operator T°[ (k) 



ff s (k) = rf s (k) + r$f(k) + ffi(k) 

= det~ 1 (A(k))[l + (g' NN - g' NA )XN(k)}Tf s (k) . (49) 

In an analogous way, the renormalized iV — > ttN decay vertex operator T^ ec (k) can be 
represented as a sum of the vacuum operator T^ ec (k), the iterated nucleon-hole T^(k) and 
A-hole T^ e A (k) contributions (see Fig. EK) - The recurrence relations for T^(fc) and Tj e A (k) 
are shown in Fig. [Ht>,c. Algebraically, they look as : 



A(k) 



rfm\ ( g'NNXN(k) 



wvj i = _ r ** (A . ) I | (50) 



Resolving (f50j) one gets the renormalized vertex operator Y^ c {k) : 



tf c (k) = Tf c (k) + r%*(k) + rfi(k) 

= det- 1 (i(fc))[l + (g' AA - g' NA ) X A (k)]Tf c (k) . (51) 

Finally, Fig. |§k shows the renormalized contact interaction —iV(k) (we drop here the 
indices of external particles for brevity) given by the sum of the bare contact interaction 
—iV(k), nucleon-hole — iVjsr(k) and A-hole — iV A (k) iterated contributions. The recurrence 
relations for —iV^(k) and —iV A (k), shown graphically in Fig. Eb,c, are : 



/l + g > NNXN(k) g' NN XN(k) 

(9'na) 



)2 XA(k) l + g' AAX A(k) 



f V N (k)\ ( g' NN x N (k)\ 

-- -V(k) 



V A (k) 



gNNXN 

g' AA XA{k) ) 



(52) 
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Resolving the system (|52|) we get a simple formula for V(k) : 



V(k) = V(k) + V N {k) + V A (k) = det" 1 (A(k))V{k) . 



(53) 



The in-medium matrix element T med is given by the sum of the direct Eq.()10|) and ex- 
change Eq. fTT| contributions where T and V are replaced by the renormalized quantities. 
The explicit formulas for calculation of |T med | 2 are given in Appendix A. 

4- Cross sections 

Fig. ITU1 shows the total in-medium pp —>■ nA ++ cross section vs the cm. kinetic energy 
above the pion production threshold 



for different densities. We compare the cross sections at fixed Q* rather than at fixed 
cm. momentum, because Q* is the direct measure of the collision inelasticity: e.g., 
a pp^tnA++iQ* = 0) = for all baryon densities, while the threshold value of the cm. mo- 
mentum is density dependent. The solid lines in the panels a,b and c of Fig. ITU1 show the 
results of the full calculation (i.e. including effective masses, the pion collectivity and vertex 
renormalization) for the combinations NLl-Setl, NL1-Set2 and NL2-Set2, respectively. We 
see that the cross section quickly drops with baryon density. The strongest effect is obtained 
for the combination NL1-Set2: in this case c^!^ n A++ drops by two orders of magnitude when 
the density increases from up to 3po- The combination NLl-Setl produces a somewhat 
weaker decrease of with density than the combination NL1-Set2. However, an 

effect of the choice of the OPEM parameters on the result is much smaller than the effect 
of the choice of the effective nucleon and A masses. The set NL2 gives less in-medium 
reduction than the set NL1 as one expects from Fig. El 

The main reason for the in-medium reduction of the cross section is the factor of 
(2m^) 3 2M^/s* in front of the matrix element squared in Eq. (}3*2*|) for the differential cross 
section. To demonstrate this, we have repeated the calculation using the NL2-Set2 com- 
bination replacing this factor by its vacuum value, i.e. (2m* N ) 3 2M^/s* — > (2m^) 3 2M A / s 
(see dashed lines in the panel c of Fig. [TUJ). One observes that the cross section calculated 
with the vacuum factor practically does not depend on the density. Further replacements of 



Q* = ys* — 2m* N — m. 



'IT 



(54) 
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the cm. momenta of the incoming and outgoing particles by their vacuum values and also 
of the effective by vacuum masses in the spectral function in Eq. have no any visible 
effect. 

In order to better understand the reason for the lack of the density dependence in the 
calculation with the vacuum phase space factor we have also calculated the cross section 
using the Set2 of the OPEM parameters with vacuum masses both in the matrix element 
and phase space factors (panels e and f of Fig. [TUJ). Panels e and f of Fig. ITU1 show the results 
with and without pion collectivity and vertex renormalization effects, respectively. Namely, 
in the calculation presented in panel f we have put all the Lindhard functions equal to zero. 
However, the A-spreading width is still present in the results shown in panels e and f. (i) As 
one can see from panel e, the cross section drops with density even without any inclusion of 
the effective masses. By comparing the solid lines in panel e with the dashed lines in panel 
c we conclude that the inclusion of the effective masses in the matrix element increases 
the cross section, (ii) From panel f we observe that the A-spreading width in the spectral 
function reduces the cross section slightly. As discussed above, this effect is comparable with 
the ambiguity given by the choice of the effective mass parameterization, (hi) By comparing 
the curves in panels e and f one sees that the pion collectivity and vertex renormalization 
reduce the cross section rather noticeably. Therefore, the density independence observed 
for the calculation with the vacuum phase space factor only (dashed lines in panel c) is, 
basically, due to a counterbalance of the effects (i) and (iii). 

Panel d of Fig. El shows the NL2-Set2 calculation with the Lindhard functions putted 
to zero in the matrix element (solid lines) and also the calculation using the matrix element 
with the in- medium nucleon and A masses from the NL2 set (dashed lines). As 
expected, both calculations are in a very close agreement since the nucleon-hole and A-hole 
Lindhard functions are absent in both cases. 

A comment is in order here on the relation to the results of Ref. where it was 

shown neglecting effective mass corrections that the in-medium cross section &™ p e ^ nA++ grows 
with baryon density. We have repeated the calculations of [15j with the parameter set B 
(A = 0.565 GeV, 

9na = 1/3) for the proton beam energy Ei ab = 0.8 GeV. Following Ref. 
I5I ] we have kept only the A-hole Lindhard function in the quasistatic approximation with 
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zero width of the A resonance 



y A (k) ~ 8f ^ k2) P " A (55) 



where wa = ffA — m^v + k 2 /(2mA). In the A spectral function the A width has been also 
put to zero jl5| . Within these approximations we have reproduced the result of Ref. 
on the in-medium enhancement of c^!^ n A++ ( see crosses in Fig. ITT| . The inclusion of the 
in-medium A width reduces the in- medium enhancement (see filled boxes in Fig. [TTf . which 
has been also demonstrated earlier in Ref. Q|. We have found that taking into account 
the nucleon-hole Lindhard function makes cr^^ nA++ a decreasing function of the density at 
p > po in the calculations with vacuum masses (filled circles in Fig. Hip. In the calculations 
with effective masses the nucleon-hole Lindhard function has a very small effect. For the 
choice of the Landau-Migdal parameters and of the cut-off factor from Set2 the cross section 
decreases faster with density (triangles in Fig. ITT]) . Introducing the effective masses amplifies 
the in-medium reduction effect strongly (see rombuses in Fig. Hip. 

Summarizing all these findings, we can state that by just changing the factor in front 
of the matrix element in Q to its in-medium value takes care of most of the in-medium 
corrections on the A-production cross section. This finding offers a way that can easily be 
implemented into any numerical realization. 

In Fig. El we show the in-medium differential cross section of the outgoing neutron at 
various densities for Q* = 0.295 GeV which corresponds to pi a b = 1.66 GeV/c at p — 0. 
For the calculation without effective mass (lower right panel) the shape of the angular dis- 
tribution does not depend, practically, on the density for 60° < O cm < 120°: only changes 
for forward/backward angles are visible. However, in the calculations with effective masses 
(all other pannels) the cross section becomes more isotropic with increasing density in qual- 
itative agreement with the Dirac-Brueckner calculations of the elastic NN cross section in 



nuclear matter 
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III. BUU CALCULATIONS 

The in-medium modified pp — > nA ++ cross section has been implemented in the BUU 
program of the version described in Ref. [27]. More details can be found in Ref. j^]. The 
cross sections for another isospin channels are related to the pp — ► nA ++ cross section by 
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the Clebsch-Gordan coefficients which gives: 

med med med med _ med ( ^,f{\ 

a pp-*pA+ ~ a pn-+pA° ~~ V^«A+ ~~ a nn-+nA° ~ g CT pp->nA++ > \ ou ) 

~med _ med 

a nn^pA- — a pp^nA++ ■ \° ' ) 

Besides the A-excitation, there are another channels of the pion production in a NN 
collision. At collision energies above the pion production threshold \fs > 2m^ + m n not 
only the A(1232), but any other resonance R can be excited : 

NN ^NR , (58) 



where R stands for N*(1440), N*(1535), N*(1650) ... (see Ref. j2J for the full list of 
resonances implemented in the BUU code). At y/s > 2{m N + m n ) the double A production 
channel opens : 

NN — >■ AA . (59) 

In addition, at collision energies slightly above the pion production threshold the s-wave 
direct pion production mechanism is important : 

NN -> NNir . (60) 

In analogy with the density modification of the A resonance cross section, we have also 
modified the processes and (joT?|) keeping, however, the constant density independent 

matrix elements as determined in Ref. 2|- 

For the in- medium cross section for the process 1)58)1 we apply Eq. (jS2J), where any higher 



resonance replaces the A resonance and \T med \ 2 is now a constant depending on the isospin 
channel. The effective mass of a higher resonance R is supposed to be 

M* R = M R + g a a , (61) 

i.e. we assume the same mass shift for all baryon resonances. The spectral function of a 
higher resonance R is given by Eq. (j2H|) with replacement "A" — ► "i?" , where the medium- 
modified total width is a sum of partial widths : 

r^(M*) = £ T R /M* R ) , (62) 

i 

rUM*) = T R /-^L , (63) 

Pi { m R) 
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where m* R = trr + g a a. m# is the pole mass of a resonance R and 

p*(M*) = J d{M B ) 2 d{M m ) 2 A B {{M B ) 2 ) 

x ^ m ((M m ) 2 ) 9(M ^f' Mm) ^(g(M-,M-,M m )//j ) . (64) 

Here, the lower indices U B" and "m" mean the outgoing baryon and meson for the i- 
th decay channel: (B,m) = {N,ir), (N,rj), (JV,w), (A, if), (A(1232), tt), (iV,p), (iV,a), 
(iV(1440),7r), (A(1232),p). A* B ((M B ) 2 ) and ^l m ((M m ) 2 ) are the spectral functions of the 
outgoing particles, Bf_ is a Blatt-Weisskopf barrier penetration factor corresponding to the 
relative angular momentum Zj of the outgoing baryon and meson. Eqs. (|62 p . ([63 )1 . (|64 )l are 
the in-medium modification of the vacuum widths from Ref. [271 ] . We have to remark 
only, that also for the outgoing A we take in Eq. (|64j) an in-medium spectral function 
AX((MX) 2 ) = <5((M^) 2 — (toa + g a o~) 2 ). This is nessary in order to keep the decay channel 
R — > Aif open in nuclear matter, since the resonance i? has now reduced mass. However, 
we did not modify any cross sections involving strange particles which is out of scope of our 
work. 

For the in-medium cross section of the process NN — > AA we apply the following for- 
mula : 

d^k^A, . (2m^ 2 (2ml) 2 



d(M 3 *) 2 d(M 4 *) 2 16ns* 



Tn 1 n 2 ^a 3 a 4 \ 2 £34 



q{Vs*,m* N ,m* N ) 



where the matrix element |T/v 1 7v 2 ^a 3 a 4 | 2 depends on the isospins of incoming and outgoing 
particles only. = 1 (1/2) if the i-th and j-th particles are in the different (same) isospin 



states. In the vacuum the cross section (J55j) is the same as in Ref. 44 1. 

Finally, we have also introduced effective mass modifications for the process of the direct 
s-wave pion production NN — > NNn which can be written in the form of the Dalitz plot 

3 = 



(2m^) 4 



A 

aa N 1 N 2 ^N 3 N 4 TT 



(27r) 3 64g(Vs*, m* N , m* N )(s*) 3 / 2 



x \T NlN2 ^ N3N4n \ 2 (y/s)S 3i d(m* u ) 2 d(ml n ) 2 , (66) 
where (m^) 2 = (pg — pi) 2 and ijn* AlT ) 2 = (pi — k) 2 are the invariant masses squared of the 



pairs of outgoing particles iV 3 , N 4 and iV 4 , tt respectively. \T NlN2 _,. N3N47 r\ 2 (^/s) is the vacuum 
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matrix element squared which depends on the vacuum cm. energy y/s (see Eq.()68|) below) 
and on the isospins of the involved particles. From Eq. one gets the following relation 
between the in-medium and vacuum cross sections : 

med / /-a g(vj» m N ,m N )s 3 / 2 fm* N \ 4 

I d(m u yd(m^ aN ^ N ^ WS) ' [b{) 

where cr^ N2 ^ N3N47T (\/s) is the vacuum cross section parameterized in Ref. 

In order to maintain detailed balance, one must modify also the inverse processes iVA — > 
NN, NR -> NN, AA -> NN, vrA^V -> NN - which are obviously essential for pion 
reabsorption - simultaneously with the direct reactions. We collect the formulas for the 
in-medium cross sections of the inverse processes in Appendix B. 

Since the mean field included in the BUU model is a Skyrme-like which is different from 
the RMF model, a recipe must be found how to match the BUU kinematics of a nucleon- 
nucleon collision with the kinematics employing Dirac effective masses. One way is to take 
the incoming particle momenta from BUU and construct s* by replacing the bare masses by 
the effective ones. However, as we discussed in Sect. II, if the cm. momentum is fixed, a NN 
collision which is subthreshold for pion production in the vacuum may get above threshold in 
the nuclear medium. This would create practical difficulties in ensuring energy conservation 
for the outgoing particles. Therefore, we have used the cm. kinetic energy above the 
pion production threshold from BUU as an input to the calculation of the in-medium cross 
sections. Namely, for each NN collision we find Q = ^ — 2m^ — m n , where \/s is calculated 
with the incoming nucleon momenta from BUU using the vacuum dispersion relation for 
nucleon energies. The in-medium cm. energy \fs* is determined by requiring that Q* = Q, 
i.e. 

\/s* = \/s — 2(m N — m* N ) . (68) 

Correspondingly, the cm. momentum of the incoming particles gets changed : 
q(y/s*,m* N ,m* N ) < q(y/s,m N ,m N ). 

Replacing only the vacuum cross sections by the in-medium ones in BUU would, however, 
not correspond to a consistent description of the resonance production/absorption rates 
because 

^12^34 OC iVEf^ (69) 
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where ^12^34 is the transition probability per unit time. 



12 - -^Q{Vs*, ml, m* 2 ) (70) 



is the relative velocity and e* = y <? 2 + ('ml) 2 , i = 1,2 are the cm. energies of colliding 
particles 1 and 2. Due to the difference between the mean fields in BUU and in the RMF 
model, the relative velocities of colliding particles are, generally, different in both models 
too. Therefore, in order to obtain the rates corresponding to the RMF model, we have 
multiplied the cross sections calculated above by the factor vl 2 /vf 2 UU , where vf 2 uu is the 
relative velocity of the colliding particles as given by BUU. 

The BUU results have been obtained by applying the soft momentum dependent mean 
field (SM, K = 220 MeV, see Refs. 13 for details) in the Hamiltonian propagation of 
baryons. Fig. shows the time dependence of the sum of the pion and A multiplicities 
in a central collision of Au+Au at 1.06 A GeV. At t = 40 fm/c practically all the A's 
have already decayed, so that the final pion multiplicity can be determined from Fig. A 
standard calculation without any in-medium effects in the cross sections ends up with 68 
pions, which is much more than the experimental value of 41 ± 7 Using the medium 
modified cross sections reduces the pion multiplicity. A bigger reduction corresponds to 
a steeper dropping effective mass with the density (c.f. Fig. EJ). Thus, the smallest pion 
multiplicity of ~ 34 is obtained within the NL1-Set2 combination. It is interesting that 
the combinations NLl-Setl and NL2-Set2 produce nearly equal pion numbers : 45 and 43, 
respectively. This is explained by the fact, that the maximum central density reached in 
central Au+Au collisions at 1 A GeV is about 2Ap , while the cross sections given by the 
NLl-Setl and NL2-Set2 combinations are similar at p < 3po (see Fig.llUj). In all calculations 
below we will always use the Set2 of the OPEM parameters, thus, we will drop for brevity 
the label "Set2". 

Fig. HU shows the average total pion multiplicity < n T > divided by the average partic- 
ipant number < A part > as a function of the total mass number A sys of colliding nuclei in 
comparison to the data [12]. The systems Ca+Ca, Ru+Ru and Au+Au were studied at the 
beam energies of 400, 1000 and 1500 A MeV. Following 12], we have fitted the total pion 
multiplicity as a function of the participant number A part by a straight line in the range 
0.15 < A part /A sys < 0.85. The participant number A part has been extracted from the geo- 
metrical overlap of colliding nuclei for a given impact parameter taking into account smooth 
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density profiles as described in [l^ . 

First, we notice that there is no a single set of parameters which is able to describe all 
data points in Fig. 1141 simultaneously. The standard calculation agrees with the data at 400 
A MeV, but it clearly overpredicts the experiment at 1000 A MeV and, for heavy systems, 
also at 1500 A MeV. 

The NL1 parameterization works for the Au+Au system at 1000 and 1500 A MeV, 
however, it is below the data at 400 A MeV. Furthermore, NL1 is unable to describe the 
experimental slopes at 1000 and 1500 A MeV. 

RHA is in the best agreement with data at 400 A MeV. Surprisingly, for the Au+Au 
system at 400 A MeV RHA gives more pions than the standard calculation. This, probably, 
can be explained by the reduced reabsorption of A's in RHA due to the in-medium reduced 
cross section AN — > NN. RHA is also rather close to the experiment at 1000 A MeV, while 
at 1500 A MeV the system mass dependence is too flat within RHA. 

The best overall description of pion multiplicities seems to be reached within the NL2 
parameterization. More precise data, in particular, at the lowest energy are needed for 
definite conclusions, however. 

Fig. 1151 shows the ratio of the average 7r~ and ir + multiplicities vs the total mass number 
of a system for Ca+Ca, Ru+Ru and Au+Au at 400, 1000 and 1500 A MeV in comparison to 
the data Q|. The calculated ratio < n w - > / < n n + > practically does not depend on the 
in-medium modifications of the cross sections, since these modifications do not change the 
ratios between the cross sections for various isospin channels. For Ca+Ca and Ru+Ru all the 
calculations agree with the data, while for the Au+Au system there is a systematic tendency 
to produce about 10 % too many 7r~'s in our calculations. Probably, the < n w - > / < n n + > 
ratio would be described better if to take into account the isospin dependence of the in- 
medium part of the nucleon propagator (jJOJ). 

Fig. shows the inclusive transverse momentum spectra of pions at midrapidity for the 
Au+Au system at 1.06 A GeV. The standard calculation (dotted line) overpredicts the data 
at p t > 150 MeV/c. NL1 (dashed line) leads to a good description at p t > 200 MeV/c, but 
it underpredicts the experiment at low p t . Almost the same result has been obtained earlier 
in using the phenomenological density dependent quenching factor for the NN <-> NA 
cross sections. A better overall description is reached within the NL2 parameterization 
(dash-dotted line), however the yields of 7r~'s and 7r°'s are still underpredicted in the region 
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of small transverse momenta. This is improved by taking into account the pion off-shellness 
in BUU [47] (solid line), which enhances the soft pion yield keeping the total pion number 
practically unchanged. 

In order to see how the density modified cross sections work in the case of light colliding 
nuclei we selected the TAPS data on 7r° production from collisions of C+C at 0.8, 1.0 and 



2.0 A GeV 



481 ] . Fig. IT7I shows the impact parameter averaged 7r° multiplicity < M >& y in a 



narrow rapidity interval near y = 0. Standard calculation is in the best agreement with the 
data for < M >Ay, while the NL1 and NL2 calculations are below the data for the beam 
energies of 0.8 and 1 A GeV. 

Finally, Fig.lTSlshows transverse mass spectra of 7r°'s at midrapidity for the C+C collisions 
at 0.8, 1 and 2 A GeV. The standard calculation provides a quite good description of the 
data [48J. In-medium modifications of the cross sections do not destroy a good description 
at 2 A GeV and also improve agreement with experiment for large transverse masses for 0.8 
and 1 A GeV, but tend to underpredict the data at small m t for 0.8 and 1 A GeV. 



IV. SUMMARY 



In this work we have calculated the in-medium NN — > iVA(1232) cross section within the 
OPEM taking into account the exchange pion collectivity, vertex corrections by the contact 
nuclear interactions and the effective masses of the nucleon and A resonance. 

(i) We have observed that even without the effective mass modifications the cross section 
decreases with the nuclear matter density at high densities if one takes into account the 
in-medium A width and includes the NN -1 Lindhard function in the calculations (Fig. 111)1 . 

(ii) Inclusion of the effective mass modifications for the nucleons and A's leads to an 
additional strong reduction of the cross section. 

(iii) The FOPI data |12| on the total pion multiplicity from the systems Ca+Ca, Ru+Ru 
and Au+Au at 0.400, 1.000 and 1.500 A GeV seem to require a dropping effective mass with 
the baryon density (NL2 version of the RMF model) in combination with the universal value 
of 0.6 for all the Landau-Migdal parameters in the spin-isospin channel. We note that this 
still depends somewhat on the size of the A-spreading width used (80 MeV p/po 0,S])- 
More precise data on the total pion multiplicity are needed in order to better select the 
model for the density dependence of the Dirac effective masses of the baryons. 
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The effect of the medium modifications of the NN <-> NA cross sections on the pion 
multiplicity depends also on the assumption about other channels of the pion produc- 
tion/absorption in NN collisions, most importantly, on the s-wave direct channel NN «-> 
NNir. Including the effective mass modifications in the NN <->• iVA channel only, does not 
lead to strong enough reduction of the observed pion multiplicity, since then more pions 
are produced in the s-wave channel. The in-medium modifications of the higher resonance 
cross sections do not influence the pion production at 1-2 A GeV collision energy sensi- 
tively : other particles like r\ and p mesons are, probably, more sensitive to higher resonance 
in-medium modifications. 

The calculations with the in-medium modified cross sections slightly underpredict the low 
transverse mass 7r° yield at midrapidity for C+C collisions at 0.8 and 1.0 A GeV (Fig. I18j) . 
This implies, that with decreasing mass number of the colliding system the in-medium 
corrections disappear faster than in our model. The reason could be a strong nonequilibrium 
momentum distribution at the beginning of collision: instead of one Fermi sphere there are 
two nonoverlapping Fermi spheres, which is not taken into account in our calculation of 
the effective mass. Indeed, the nonequilibrium momentum distribution reduces the scalar 
density, and, therefore, increases the effective mass. This effect is expected to be stronger 
for lighter systems, since smaller baryon density is reached in this case. 

Another open problem is the general tendency of the transport models to overpredict 
pion multiplicity in heavy colliding systems at AGS energies (c.f. Ref. 111]). This problem 
can not be solved by just density modifications of the resonance production/absorption 
cross sections, since most pions are produced by the string mechanism here. In-medium 
modifications of the FRITIOF model, thus, could help to solve this problem. 



Acknowledgments 

Numerous stimulating discussions with M. Post and Dr. L. Alvarez- Ruso are gratefully 

acknowledged. The authors are grateful to Markus Post for the reading of the manuscript 

before publication, critical comments and suggestions and to Luis Alvarez-Ruso for pointing 
I I 

out Ref. 28j which has influenced our results. 



22 



APPENDIX A: CALCULATION OF \T med \ 2 FOR THE PROCESS pp -> raA +H 



For the calculation of |T med | 2 it is convenient to introduce the renormalized coupling 
constants 

f A (k) ee det- 1 (i(A:))[l + (g' NN - ^ a )xat(*)]/a(* 2 ) , (Al) 
f {k) = d et- 1 (i(A;))[l + (g' AA - g' NA ) X A(k)]f(k 2 ) (A2) 

and the renormalized Landau-Migdal parameter 

a( k ) = det(A(fc))^ A 

9{ ' ~ [1 + Wnn - g'NA)XN(k)] [1 + (g' AA - g' NA ) X A(k)] ' 1 J 

In terms of f\{k), f(k) and the renormalized vertex operators and the renormalized 

contact interaction have the same form as the vacuum ones : 

ff s (k) = -^®(s . k)T 7 , (A4) 
tf\k) = ^(cr-k)^ , (A5) 
= fA{k) { {k h (k)[A(^S a Tx(2)}[x\S)a a r X (l)} . (A6) 

Notice, that /a(^) 5 /(&) and g(/c) are complex-valued quantities. Substituting (|A4j) . ()A5|) 
and ()A6|) into (fTUj) on the place of corresponding vacuum quantities and also replacing 
(S • k) — >• (S ■ k 4 ) (see Eq. ()21|l ) in Eq. (jA4|) one can get after some algebra : 



1 l r rmed]2 _4\fA(k)\V(k)\ 



x k 4 V|D(A ; )| 2 + (r(A:)D(A ; ) + ^(^*(fc))(k 4 .k) + 3|^(A;)| 2 ] , (A7) 

where the term — t/k 2 is introduced in order to restore the Lorentz invariance, as discussed 
in Sect. II. The exchange term \T™ ed \ 2 j s given by Eq. (|A7|) with replacements k — > k', t — > u 
and k 4 -> k^ (see Eq.(|22|). 

For the spin-averaged interference term we get : 

f*A(k')r(k')f A (k)f(k) 

Ai,A2,A3,Aa4 '^ 7r 



X 



{l[2(k^ • k 4 )(k' ■ k) + (k 4 ■ k)(ki ■ k') - (k 4 • k')(k' 4 • k)]D*(k')D(k) 

(k 2 



+ ~(k 4 ■ k')D*(k')g(k) + ^(k 4 • k)D(k)~g*(k') + 4? WW)} + cc. . (AS) 
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APPENDIX B: IN-MEDIUM CROSS SECTIONS FOR THE INVERSE PRO- 
CESSES 



We use below the detailed balance principle, i.e. the equality of the matrix elements of 
the direct and inverse processes, in order to get the cross sections for the reactions iVA — > 
NN, NR — > iViV, AA — ► iViV and the pion s-wave absorption rate by two-nucleon pairs 
itNN -> NN. 

The in-medium differential cross section for the process iVA — ► NN is given by the 
expression (c.f. Eq. (J32J)) : 

d^: d A^ Nl N 2 (V^) = (2m* N f2M* A 2 ^- ^(y^, m* N , m* N ) 

dSl 64ttV 2J a + 1 1 q(y/s^,M A ,m* N ) 12 ' 1 ' 



where \T med \ 2 is calculated in Appendix A (Eqs. (|A7[) . ([A8p ). J a = 3/2 is the spin of the 
A-resonance. For the process NR —+ NN, where R means one of the higher resonances, 
one should simply replace in Eq. (|Bl|) "A" by "i?" : now and Jr are the effective mass 



and spin of the resonance R, respectively. \T med \ 2 is now a constant depending only on the 
isospin channel. 

The in-medium cross section of the reaction A A — > NN reads as follows (c.f. Eq. (J65|) : 

(2m* N ) 2 (2m* A ) 2 4 



med — \ N) \*""'AJ IT 1 2 



X ^^ m N^N) s (B2) 



The pion s-wave absorption rate by the two-nucleon pairs is : 

™ e(!>/ n 1 f 9d 3 Pi _ gd 3 p 2 gd 3 p 3 M _ , gd 3 p 



{k) = 2T ^J (2^ ni (2^ n2 (2^ 
x (2^5^ (k + p{ +p* 2 -pi -pl)S 12 S u (2m* N y\T N3N4 _ NlN2n \ 2 ( v G) , (B3) 

where g = 2 is the nucleon spin degeneracy factor, n,i = n(pi), (i = 1,2,3,4) are the 
nucleon occupation numbers. The sum in Eq. (jB3|) is taken over the isospin projections of 
the incoming nucleons 1 and 2. The vacuum matrix element squared 

TTf Hi rz\ (27r) 3 64g(y^, m N , m N )s^ 2 a v ^ Ni ^ NlN2 ^) 

\1 NaNi-^NiNw \W S ) = 777 7T75 TTl \2M \2 ^ 

depends on the vacuum cm. energy and on the particle isospins. The integral ()B3)1 is 
performed by the Monte-Carlo method as described in Ref. j^J. 
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Table 1. Parameters of the medium modified OPEM 



Set # 


9nn 


9na 


9aa 


A, GeV 


f 


0.9 


0.4 


0.4 


0.70 


2 


0.6 


0.6 


0.6 


0.67 
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FIG. 1: Direct (a) and exchange (b) diagrams contributing to the amplitude of the process 
N1N2 —>■ N3A4. The wiggly line denotes either ir exchange or the contact interaction oc g^A- 
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FIG. 2: Upper panel: the total cross section pp — ► nA ++ in vacuum as a function of the cm. 
energy y/s. Calculations are performed using the Set 1 (dashed line) and Set 2 (solid line) of the 
OPEM parameters (see Table 1). Results obtained within the model of Ref. ly] are shown by 
dotted line. Lower panel: vacuum cross section of a single 7r + production in pp collisions vs y/s. 
. Results with the pp — * nA ++ cross section calculated with the Set 2 and using the model of Ref. 
1^ | are shown by solid and dotted lines respectively. The s-wave background, isospin 1/2 and 
3/2 (excluding ^33(1232)) resonance contributions \l\ are shown by short-dashed, long-dashed and 
dash-dotted lines respectively. Solid and dotted lines show the incoherent sum of the pp — ► nA ++ 
cross section multiplied by 10/9 and the background plus isospin 1/2 and 3/2 higher resonance 
contributions. The data are from Ref. 




FIG. 3: Outgoing neutron cm. polar angle distribution in the reaction pp — * npir + at pi a i, = 1.66 
GeV/c. In the calculations only the pp — ► nA ++ channel is taken into account. The curves are 
labeled as in the upper panel of Fig. [2 The histogram shows the data from Ref. 
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FIG. 4: Invariant mass distribution of outgoing tt + and proton for pp collisions at pi a b = 1.66 
GeV/c. In the calculations only the pp — ► nA ++ , A ++ — ► 7r + p reaction channel is taken into 
account. The data from Ref. 30] are shown by the histogram. 




FIG. 5: Ratio of the effective nucleoli mass to the vacuum nucleon mass vs the ratio of the baryon 
density to the nuclear saturation density Dashed, solid and dash-dotted lines show the calculations 
using NL1, NL2 and RHA respectively. 




FIG. 6: (a) - Pion self-energy —ill(k). (b),(c) - The recurrence relations for the AiV -1 and 
AW -1 iterated contributions (see Eqs.(|H4|)). Arrows are not shown on the AiV -1 loops, since 
each of these loops is the sum s- and it-channel graphs (with anticlockwise and clockwise directed 
arrows, respectively). 




FIG. 7: (a) - Renormalized irNA vertex operator f " (k). (b),(c) - Recurrence relations for the 
AA^ 1 and NN^ 1 iterated contributions to ff s (k) (see Eas.pg))). 




FIG. 8: (a) - Renormalized nNN vertex operator Tj(k). (b),(c) - Recurrence relations for the 
AA^ 1 and NN^ 1 iterated contributions to tf c (k) (see Eas.(|5D|l). 




FIG. 9: (a) - Renormalized contact interaction —iV{k). (b),(c) - Recurrence relations for the 
AiV -1 and NN^ 1 iterated contributions to —iV(k) (see Eas.(|52j0. 
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FIG. 10: In-medium pp — > nA ++ cross section as a function of the kinetic energy above the pion 
production threshold Q* (see Eq.()54[)) for various parameter sets (see text for details). In all cases, 
excepting the dashed lines in panel (c), lines from the uppermost to the lowermost correspond to 
the baryon densities p = 0, po, 2po and 3po, respectively. For the dashed lines in panel (c) the 
density order is opposite. 
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FIG. 11: Density dependence of the ratio of the in-medium and vacuum total cross sections 
pp — > nA ++ for Q* = 0.225 GeV (in vacuum this corresponds to the proton beam energy of 0.8 
GeV). Line marked with crosses corresponds to the approximations of Ref. 2] with the parameter 
set B [l5j. Boxes and circles show the calculations including the in- medium A width without and 
with the NN^ 1 Lindhard function contribution respectively. Triangles and rombuses show the 
results for the Set2 of the OPEM parameters without and with inclusion of the effective mass, 
respectively. 




FIG. 12: Neutron cm. polar angle dependence of the in-medium differential pp — > nA ++ cross 
section. The lines from the uppermost to the lowermost correspond to the baryon density p = 
0, po, 2po and 3po, respectively. Calculations are performed for Q* = 0.295 GeV. 




FIG. 13: Time dependence of the sum of the pion and A multiplicities in the central collision 
of Au+Au at 1.06 AGeV. The standard BUU calculation is shown by the dotted line. The dash- 
dotted, solid, short-dashed and long-dashed lines show the results with the in-medium modified 
resonance production/absorption cross sections by using the combinations RHA-Set2, NL2-Set2, 
NLl-Setl and NL1-Set2, respectively. 
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FIG. 14: Total pion yield per participating nucleoli for different energies and colliding systems. 
Dotted, dash-dotted, solid and dashed lines correspond to the standard, RHA, NL2 and NL1 
calculations, respectively. The data points are from Ref. |12|]. The results of the data analysis A 
(B) are represented by filled circles (stars) 




FIG. 15: Ratio of the charged pion yields for different energies and colliding systems. Various 
curves are denoted as in Fig. 1141 The data points are from Ref. 
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FIG. 16: Pion inclusive transverse momentum spectra at midrapidity for the reaction Au+Au 
at 1.06 AGeV. Dotted lines show the standard BUU calculation. Dashed and dash-dotted lines 
show results with in-medium modified resonance production/absorption cross sections applying the 
NL1 and NL2 parameterizations, respectively. Solid lines correspond to the NL2 parameterization 
taking into account pion off-shellness. The normalized rapidity intervals in which spectra are 
extracted are Y<® = —0.2 -j- 0.2 for vr* and = -0.25 0.21 for tt°, where = (y /y pro j)c.m, 
The data are from Ref. Il2|. 




FIG. 17: The average tt multiplicity in a narrow rapidity interval near midrapidity as a function 
of the beam energy for C+C collisions. Dotted, dashed and solid lines correspond to the standard, 
NL1 and NL2 calculations, respectively. The rapidity intervals are y = 0.42 -j- 0.74 at = 0.8 
and 1.0 A GeV, and y = 0.80 -=- 1.08 at E lab = 2 A GeV. The data are from Ref. |4j|. 
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FIG. 18: Transverse mass spectra of tt° mesons in the narrow rapidity intervals near midrapidity 
(see caption to Fig.HZJ for the 12 C+ifC system at 0.8, 1.0 and 2.0 A GeV. Curves are denoted as 
in Fig. 1161 The data are from Ref. 



